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Abstract 

Radial solitons are investigated in armchair carbon nanotubes using 
a generalized Lennard- Jones potential. The radial solitons are found in 
terms of moving kink defects whose velocity obeys a dispersion relation. 
Effects of lattice discreteness on the shape of kink defects are examined 
by estimating the Peierls stress. Results suggest that the typical size for 
an unpinned kink phase is of the order of a lattice spacing. 

1 Introduction 

Nanostructured materials are known for their unique properties characterized 
by rich mechanical, structural and electronic properties on a short length 

scale. These particular properties give promizing future for various technolog- 
ical applications, and have motivated great recent interest to those materials. 
Among the nanostructured materials recently investigated, carbon nanotubes 0] 
appear as the most fascinating due in part to their cylindrical structure CI] ■ 
Raman spectroscopy El has been quite helpful for the understanding of their 
structural properties, e.g. by allowing direct measurements of their vibration 
spectra. On the one hand, experimental results show that the interplay of 
the structure and the short length scale gives rise to a marked dependence of 
the frequency of some vibration modes on the tube diameter as for instance 
the so-called A g breathing mode. On the other hand, combined Raman spec- 
troscopy and inelastic neutron scattering experiments reveal pronounced elastic 
properties as well as plastic deformation processes, two foundamental features 
suggesting that nanotube involve many-body interactions. Owing to their weak 
and dispersive characters, many-body interactions are most often non-local or 
non-bond forces with van der Waals features IB] ■ Therefore they can readily 
be thought of as the main source of stability of nanotubes as well as nanotube 
bundles such as fullerenes. 

A large amount of work has recently been devoted to the vibration spectra of 
nanotubes. Being radial, the breathing mode emerges as the most influenced by 
the tube geometry 6, 8, 9 . This feature is manifest in its vibration frequency 
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which displays a strong dependence on the tube diameter. Theoretical attempts 
to explain this dependence was carried out following pair-potential approaches 
and led to a rather simple relationship [5] El agreeing with experiments. More 
explicitely, in these theories many-body potentials are used to deal with the 
dispersive interactions. The vibrational spectrum is thus constructed in terms 
of the spectrum of phonon excitations in the presence of the dispersive interac- 
tions. In general, two classes of phonons show up in the Raman spectra namely 
tangential and radial phonons. 

The recent interest to excitations in carbon nanotubes was also directed to- 
ward modes involving amplitudes more larger than phonons. Solitons, which 
are robust objects describing these large-amplitude excitations, have so far been 
suggested in two main distinct ways El El The first El relates to the sponta- 
neous lattice distorsion responsible for the dimerization of the electronic ground- 
state of the nanotube. In fact, this first kind of soliton excitation is well known 
in systems undergoing Pcicrls instabilities, where it forms the state of broken 
symmetries then called domain wall. The second kind of soliton El follows from 
the Hamiltonian description of the nanotube dynamics, i.e. by writing down a 
total Hamiltonian consisting of a kinetic energy which describes the tube dy- 
namics, and a potential energy accounting for the many-body interaction. The 
potential energy in this last case connects to the Brenner potential El; which 
in the particular context of the paper is expanded for weak interatomic 
displacements. Interestingly, this expansion leads to a generalized nonlinear dif- 
ferential equation which is just the so-called nonlinear Klein-Gordon equation, 
and which can in principle be transformed into any of the known dispersive 
nonlinear equations as the Nonlinear Schrodinger equation, Korteweg-de Vrie 
equation El > etc. Soliton solutions of these dispersive equations are particularly 
interesting for the dynamics of atomic lattices because of their robustness due 
to the balance of dispersion by nonlinearity both coming from the many-body 
potential. 

In this work, we will attempt to formulate the Hamiltonian description of the 
soliton excitation in single- wall nanotubes in terms of radial solitons. To remain 
as close as possible to the theory developed in Eh we keep the spirit of the 
many-body interaction generating both the dispersion and nonlinearity in the 
system dynamics. However, unlike this work we will use a generalized Lennard- 
Jones E1E] potential hereafter called "Lennard- Jones hierarchy" (LJh). Let us 
remark that a simple thoeretical model for such an Hamiltonian description fol- 
lows from the assumption that the carbon-carbon interaction, readily regarded 
as the source of the many-body effects, promotes a ID atomic lattice structure 
in which atomic sites interact via the many-body potential. In single-wall nan- 
otubes, several of these II? atomic lattices result into a nanotube sheet and the 
effect of some rigid background substrate that may provide nonlinearity to the 
system dynamics as in other contexts El E] is rather unlike. In fact, being 
a nonlinear function of the lattice displacement, the many-body potential will 
furnish the necessary dispersion and nonlinearity needed to stabilize the lattice 
excitations. 

Competing effects between dispersion, nonlinearity and lattice discreteness can 
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become quite relevant to the existence and stability of some of these nonlinear 
excitations. For instance, kink-soliton defects will show a rather strong sensibil- 
ity vis-a-vis the discrete lattice structure and will most often need a threshold 
width to survive lattice discreteness effects. To this viewpoint, it is well es- 
tablished that narrow kinks are more prone to form pinned soliton states |2(J| . 
proving that pinning effects can have dramatic incidence on transport properties 
of materials |21U22ll2"3*| . As we will be interested with kink-soliton defects, it will 
be instructive to look at the effects of the comptctition between the two main 
ingredients of the model i.e. the dispersion and nonlinear, on the continuum 
shape of a kink-soliton defect propagating on the discrete ID lattices forming 
the single- wall nanotube. 

To start, we investigate the dynamics of radial kink-soliton defects in the arm- 
chair nanotube for many-body potentials belonging to the LJh. We will also 
examine dispersion features of the continuum kink-soliton solution in the ID 
discrete lattice and derive a dispersion law, valid for a given (m, m) armchair 
configuration. Next, the effects of lattice discreteness on shapes of the radial 
kink-soliton defects will be considered. In this context, the analytical expression 
of the Peierls stress will be derived. 



2 The model, soliton solutions and dispersions. 

Consider a ID lattice consisting of an atomic backbone in which neihgbouring 
atoms interact via a pair potential of one of the LJhs (n\, 712). Denoting by r n 
the relative displacement of the n th atom, the total energy of the lattice can be 
written: 

M 

— f% + V LJh (r n+1 - r n ) 



(1) 



where dot refers to derivative with respect to time t, M is the total mass asso- 
ciate to the n th site and Vljk is the LJh potential assumed of a general form: 



V L jh{r) =4e c [(^) -2(7 



a \n 2 



(2) 



ni and ni in this last formula are integers chosen such that n\ > 112, while a 
is a characteristic length and e G a characteristic energy. The LJh potential J5J 
possesses two different interaction branches i.e. a short-range repulsive branch 
and a long-range (van der Waals-like) attractive branch. The two branches meet 
at the equilibrium point: 



<t, V L jh{r ) = -4 — 
2n 2 



(3) 



corresponding to the dissociation point of the lattice. For the bare LJ potential, 
(ri2,ni) = (6,12) suggesting a potential minimum at r = a, for which V eq = 
— 4e . Another frequent set of values |17j is {n2,n\) = (3,9) for which r Q = 



(|) 1/6 a and V et 



.1/2 
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Figure 1: Sketch of the one-dimensional lattice of atoms in an armchair nan- 
otube. u is a unit vector. 

The characteristic distance a will generally depend on geometry of the material. 
Nevertheless, by appropriate means we can also define a square frame in an 
hexagonal geometry and parametrize the displacement vectors r n by defining 
radial variables 9 n with respect to the new frame, obtaining: 

r n+ i - r„ = a n +iu n +i - aiU n , (4) 

where u n are unit vectors along rods lying on the sheet of nanotube as in figure^ 
while resting on atoms(labeled i, i + 1, ...) of the substrate backbone. At 
equilibrium, the rods align uniformly with an equilibrium angle 9 a with respect 
to the ID lattice. The equilibrium angle 9 can thus be assumed as the chiral 
angle of the nanotube, and for an armchair configuration 9 a will be determined 
by the "chiral coordinate" (m,m) which also fixes the diameter of the tube. 
Therefore, we readily expect this last quantity to affect excitations of the tube. 
A simple way to implement this dependence is to consider instead the twist 
motion of the entire rod then assumed rigid. In this way, the mass M turns to an 
effective mass resting on a reference atom on the substrate lattice. Expressing 
Uji as a function of the radial coordinates and discarding the constant chiral 
angle, we obtain: 

u n = e x cos 6 n + e y sin 9 n . (5) 

Since we assume rigid rods, we can set a n +i = a n — ^o- Furthermore, we 
place ourselves on an orthonormal basis where e.i =x ^ v are orthonormalized unit 
vectors. From (JjjJ and taking into account all the considerations made above, 
the separation between the (n + l) th and the n th site block along the ID lattice 
reduces to: 

r 2 = 2^[l-cos(0 n+ i-0„)] (6) 
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We will now proceed to the key transformation of our approach, consisting of 
an expansion of in Taylor series arriving at: 



k=l 



V L jh = V (n 1 ,n 2 ) + ^ Vfc(wi,n 2 ) cos fc [6 n+1 - 6 n ) 

where we have set: 

V k (ni,n 2 ) = 



•>i) ( a 



2b 



(n 2 ) / cr 



uU) 



k \V2£ J k \V2i c 

(-l) fe j(j + 2) (j+2k-2) 



(7) 



(8) 



(9) 



* fc! 2 2 2 
This expansion is a simple generalization of the Taylor series representation 
of the LJ potential, except that the condition of weak relative displacements 

' n+l 'n 

is replaced by its equivalent involving the relative angle 6 n +i ~ n . 
Similar expansion has also been used in the study of critical properties of lipid 
monolayers |27| . 

When fc = 1, equations (7J- lead to a one-period sinusoidal potential i.e. the 
sine-lattice(SL) potential |281 1291 1501 1311 132| . To complete our transformation, 
we also need to formulate the kinetic energy in terms of the discrete radial 
variables. This is achieved by introducing the elementary displacement dr n ~ 
£ d9 n . With this assumption, the total Hamiltonian Q reduces to the SL 
Hamiltonian i.e.: 



e S l = £ -^el + v SL (0n+i - e n ) 

n=l L 

Vsl{6u+i - Qn)} = Vi{n 1 ,n 2 ) cos(0 n+ i - 9 n ) 
from which we derive the following set of difference differential equations: 



with: 



i+i - n ) + sin (■ 
2 ^1(^1,^2) 



l -x-e n )] = o 



Mil 



(10) 



(11) 
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Equation i|12|) . also called SL equation 28, 29, 3Tni3"Tll32| . was first introduced for 
conformational transitions in DNA macromolecules. It admits, at large radial 
displacements, a kink-soliton solution described by the ansatz: 



n (t) — 2 arctanexp (qna — ujt) 



(13) 



where a is the equilibrium separation between two neighboring atomic sites 
along the ID lattice. The ansatz (|13H can also be looked on as a kink-soliton 
defect that propagates along the ID lattice being subjected to the dispersion 
law: 



uf = Auo 2 sinh (qa a /2) 



(14) 
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In this last relation, the quantity lo q acquires more physical meaning in terms 
of the characteristic frequency of kink vibrations on the ID lattice. According 
to the argument of the characteristic width of the kink defect is inversely 
proportional to q while the kink velocity is determined by the dispersion rela- 
tion Q14p. To this viewpoint, from 114|l a group velocity can indeed be defined 
as follow: 

This group velocity displays a threshold value in the limit q —> and is given by 
-d = ui a . Returning to formula (|12|) . we will notice that ignoring the constant 
Vi (721,712) which depends mainly on details of the many-body potential, the 
characteristic frequency oj appears as an inverse function of the tube diameter 
l a . It turns out that the threshold group velocity # will dependent on the ratio 
between the two characteristic length scales of the nanotube i.e. the lattice 
constant a and the tube diameter £ a . While a can be assumed a free parameter, 
to is a function of the chiral coordinate (to, to) meaning that the ratio a /£ will 
acquire distinct values for different armchair configurations. 
On figure |5J we plot the reduced group velocity d{q)/ n o as a function of the 
kink width 1/q. We observe that the kink velocity gets slowed down when its 
width increases. Moreover, the group velocity rapidly saturates to the minimum 
threshold -d at finite, but relatively large values of the kink width. It is useful to 
underline the qualitative agreement between the behaviour reflected by curve of 
figure|3and result of ref. ^3] on longitudinal solitons(see figure 2 in this paper) 
which we recall it, uses the Brenner potential. 

We will close this section by looking at the energetics of the radial soliton defect 
obtained in (|f 3(1 . In particular we estimate the defect energy which in the usual 
language, is the characteristic energy needed to create a kink soliton of this 
form. Replacing in the continuum version of the SL Hamiltonian i|l(J|) ('i.e. 
making the continuum limit approximation na Q — > x) and taking into account 
the dispersion relation l(14j) . we find: 

E=^- (16) 
qa Q 

This quantity is vanishing at large defect widths, a behaviour already displayed 
by the kink frequency or also the dispersion relation (|I4|) . 



3 Discrete lattice effects on the kink-soliton de- 
fects. 

The analytical expression (|I3[) is actually not an exact solution of the nonlinear 
equation but describes a continuum soliton solution accomodating the ID 
discrete lattice. Then, in addition to the dispersion relation l|14[) which dictates 
the propagation velocity of the continuum kink soliton, the discreteness of the 
propagation medium will introduce an additional constraint for the stability of 
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Figure 2: The reduced group velocity of the kink defect as a function of the 
dimensionless soliton width(i.e. in unit of the lattice constant a Q ). 
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shape of this continuum kink soliton in the II? discrete lattice. This constraint 
acts as a mechanical stress that tends to keep the soliton pinned to the discrete 
background of the lattice, hence it is comparable to the Peierls stress commonly 
observed in crystal monolayer as well as multilayer phenomena |33) . The concept 
of second Peierls stress is best known for systems involving a rigid background 
substrate, such as in the Frenkel-Kontorova JS| ancl Frank-van der Merwe 
models, in addition to the strain energy of a Hookes type. Nevertheless, in a 
recent paper [2S| we suggested an approach to deal with the Peierls-Nabarro 
effect for model Hamiltonians of the SL type and will follow this work. Start- 
ing, remark that taken in its current analytical form i.e. keeping only the site 
coordinate n, the continuum kink soliton (|13|) looks strongly dependent on the 
discreteness and will not survive pinning to the lattice sites. To avoid this 
dramatic lattice discreteness effect, we shift the kink position by introducing a 
center-of-mass coordinate n that preserves both the kink translational-invariant 
shape and dispersion governed by the dispersion rclation lfHEjl . With these con- 
siderations the kink shape now becomes: 

9 n (t) — 2 arctanexp [q (n — n )a ] (17) 

Inserting l|17l) in the strain term of the Hamiltonian I|1U|) and Fourier transform- 
ing as in [26], we obtain: 

p a q 2 sinh(2n 2 / qa ) 

e p , or the Peierls-Nabarro potential barrier for the SL model l|10ll . is plotted on 
figure versus the soliton width. Here we use the quantity: 

V = 87T 2 Ma o e 2 uj 2 . (19) 

as unit of mechanical stress. The physical meaning of the decrease of the Peierls- 
nabarro potential barrier with an increase of the soliton width is a reduced effect 
of the mechanical stress for relatively large soliton defects. Typically, a soliton 
defect of size about one or a few lattice spacings will almost be insensitive to 
the lattice discreteness as suggested by the drastic fall-off in figure [3] What this 
result implies as for the physics of the nanotube is the possibility of a conducting 
state dominated by narrow soliton defects. 



4 Summary and concluding remarks 

We have explored an Hamiltonian description for radial soliton excitations and 
dispersions in single-walled nanotubes with an armchair configuration using a 
generalized Lennard- Jones potential. A first point to underline is the agree- 
ment between our results and a previous study |13| on longitudinal solitons us- 
ing Brenner's potential. This includes the velocity and energy whose variations 
with the soliton width are qualitatively similar to predictions of this previous 
work. In addition, we found that dispersion confines the soliton velocity within 
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Figure 3: Amplitude of the Peierls stress, plotted versus the dimensionless soli- 
ton width(i.e. in unit of the lattice constant a ). 
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a defined interval characterized by a minimum threshold velocity proportional 
to the ratio of the lattice spacing to the tube diameter. 

Effects of the lattice discreteness on the soliton shape have also been investi- 
gated. We obtained that the threshold width for an unpinned soliton phase was 
nearly of the order of a lattice spacing. This narrow soliton-induced unpinned 
kink defect phase is the signature of a dominant dispersion and may therefore 
well be understood as reflecting the absence of a single-particle potential. In- 
deed, in general the single-particle potential will enhance the threshold width 
for the unpinned kink phase [25] owing to the rigid character of the substrate to 
which it relates. Whether this single-particle substrate potential is pertinent to 
single- wall nanotubes and/or their derivatives i.e. nanotube bundles, fullerenes, 
etc., is to our view a central issue for a more realistic description of the nonlinear 
dynamics of these systems. 
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